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I.  INTRODUCTION 

1.  General.  This  report  describes  a technique  for  measuring  energy  losses  in 
type  11  superconductors  due  to  alternating  current  iiow  in  superconducting  solenoids. 
Superconductivity  and  the  transient  loss  theory  are  briefly  discussed,  and  the  apparatus 
used  to  measure  the  losses  is  described  in  detail.  Loss  data  as  a function  of  frequency, 
which  was  obtained  from  measurements  on  a number  of  solenoids  wound  from  com- 
mercially available  multifilament  (MF)  NbTi  conductors,  is  presented. 

2.  Background.  A class  of  metallic  conductor  exists  which  when  cooled  to  very 
low  temperatures  loses  all  electrical  resistance.  Such  metals  are  called  superconductors. 
This  phenomena  was  first  discovered  in  1911  and  was  accompanied  by  glowing  predic- 
tions of  widespread  application;  however,  it  was  soon  discovered  that  this  resistance- 
less  property  was  destroyed  by  very  small  magnetic  Fields,  eliminating  typical  electro- 
magnetic applications. 

In  I960,  high-field,  type  II  superconductors  were  discovered  by  Kunzler. 1 
These  superconductors  remained  in  the  superconducting  state  when  exposed  to  high 
magnetic  fields  and  high  currents,  and  again  superconductivity  was  considered  as  hav- 
ing a bright  future  of  industrial  and  military  applications.  It  is  these  type  II  supercon- 
ductors which  are  the  subject  of  the  research  described  in  this  report.  In  spite  of  the 
favorable  electromagnetic  properties  exhibited  by  these  materials,  no  widespread,  prac- 
tical applications  for  them  have  been  found  in  the  16  years  since  their  discovery. 

3.  Superconductivity.  The  resistance  of  materials  in  the  superconducting  state 
is  zero  for  all  practical  purposes.  Recent  measurements  have  placed  an  upper  limit  of 
10'21  ft  on  the  resistivity  of  a superconductor.  This  compares  with  a value  of  10'9  ft 
for  high-purity  copper  at  4.2K.  The  critical  temperature,  Tc,  is  the  temperature  at 
which  a given  material  undergoes  the  transition  from  the  normal  to  the  superconduct- 
ing state,  but  it  is  the  magnetic  behavior  of  superconductors  which  divides  them  into 
classes  known  as  type  I or  type  II. 

A bulk  specimen  of  metal  in  the  superconducting  state  exhibits  perfect  dia- 
magnetism with  the  magnetic  induction  B = 0.  This  is  called  the  Meisner  effect.  In 
reality,  a magnetic  field  can  penetrate  a small  distance  into  a bulk  superconductor.  In 
the  pure  superconducting  state,  the  only  field  allowed  is  exponentially  damped  as  it 
goes  in  from  an  externa)  surface  with  the  strength  of  the  field  decreasing  inward  from 
the  surface  according  to  the  relation 

* J.  E.  Kunzler,  Rev.  of  Mod.  Phy.,  Vol.  33,  No.  1,  1961. 
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PI  T.  ^ 


B(X)  = B(0)e  - X 
X 


where  X is  the  distance  in  from  the  surface  and  X is  a characteristic  length  known  as 
the  penetration  depth.  This  distance  is  quite  small  a typical  value  being  5.0  x 10'6 
cm  for  tin.2 

A sufficiently  strong  magnetic  field  will  destroy  superconductivity.  This 
threshold  on  critical  value  of  the  applied  magnetic  field  is  denoted  by  Hc  and  is 
referred  to  as  the  critical  field.  This  value  is  dependent  on  temperature,  but  the  value 
commonly  used  is  either  T = 0 for  pure  superconductors  or  T = 4.2  for  practical  super- 
conductors. 

There  is  an  upper  limit  to  the  current  which  can  be  passed  through  a super- 
conductor without  driving  the  material  back  into  its  normal  condition.  There  is  a 
close  relation  between  this  critical  current,  lc,  and  the  magnetic  properties  of  the 
superconductor. 

4.  Type  II  Superconductivity,  There  is  no  difference  in  the  fundamental 
mechanism  of  superconductivity  in  type  I and  type  II  superconductors;  however,  the 
Meisner  effect  is  quite  different  in  the  two  types.3  A type  I superconductor  (Figure  1 ) 
excludes  a magnetic  field  until  superconductivity  is  destroyed  suddenly  and  com- 
pletely, and  then  the  field  penetrates  completely.  A type  II  superconductor  (Figure  1) 
excludes  the  field  completely  only  in  relatively  weak  fields  up  to  a field  Hc,  called  the 
lower  critical  field.  Above  Hc,  the  field  is  partially  excluded,  but  the  material  remains 
electrically  superconducting.  At  a much  higher  field,  typically  6T  to  120T,  the  flux 
penetrates  completely  and  all  superconductivity  vanishes.  This  field  is  called  Hc2  - 
the  upper,  critical  field.  It  is  this  ability  of  type  II  superconductors  to  remain  elec- 
trically superconducting  at  high  magnetic  fields  that  makes  them  of  considerable 
interest  for  practical  applications. 

In  the  region  between  Hcl  and  Hc2,  a type  II  superconductor  is  said  to  be 
in  the  vortex,  or  mixed,  state.  In  this  state,  flux  penetrates  the  specimen  as  an  array 
of  flux  tubes  parallel  to  the  field.  Each  flux  tube  consists  of  a normal  core  surrounded 
by  a screening-current  vortex.  As  the  field  is  increased,  the  number  of  flux  tubes  in^ 
creases  until  at  the  upper,  critical  field  the  normal  regions  overlap,  and  supercon- 
ductivity is  destroyed.  The  flux  motion  leads  to  hysteresis  in  the  magnetization  of 
type  II  materials  and  is  a major  source  of  power  losses  under  alternating  current 
conditions. 


2 _ , , _ 

G.  G.  Haseldenj,  ed..  Cryogenic  Fundamentals,  Academic  Press,  1971. 

C.  Kittel,  Introduction  to  Solid  State  Physics,  4th  ed.,  John  Wiley  & Sons,  1971. 
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Most  useful  superconductors  for  technological  purposes  are  type  II  materials. 
Superconductors  in  this  class  have  the  highest  values  of  critical  temperature,  field,  and 
current  and  thus  offer  considerable  technological  advantage  over  type  I materials. 
Throughout  the  remainder  of  this  report,  the  word  “superconductor”  will  be  used  to 
mean  type  II  superconductor  unless  otherwise  specified. 

5.  Practical  Superconductors.  Currently,  there  are  several  alloys  used  to  make 
commercial  superconductors.  The  most  commonly  used  are  niobium  titanium  (NbTi) 
and  niobium  tin  (Nb3Sn)  which  are  integrated  in  a matrix  of  good,  normal  conductor 
usually  copper.  This  matrix  is  essential  to  produce  a stable,  current-carrying  wire.4  5 

Niobium  titanium,  which  has  a critical  temperature  of  10K  and  an  upper 
critical  field  of  1 2T,  has  the  advantage  of  being  very  ductile  so  that  fabrication  is  easier 
and  the  superconductor  cheaper.  Reducing  the  size  of  the  current-carrying,  super- 
conducting filaments  reduces  hysteresis  losses  and  leads  to  more  stable  characteristics 
under  transient  field  conditions.  Because  of  its  high  ductility,  NbTi  can  be  easily 
manufactured  in  MF  conductors.  These  are  wires  with  many  small  superconducting 
filaments  (Figure  2).  The  NbTi  filaments  have  diameters  typically  in  the  range  of 
1-30  micrometers  (/am),  and  over  1000  filaments  can  be  incorporated  in  one  wire 
approximately  0.015  inch  in  diameter.  These  filaments  are  usually  twisted  inside  the 
matrix  so  that  they  spiral  along  the  length  of  the  conductor  with  a twist  pitch  of  1-10 
per  inch. 

Twisting  is  especially  important  in  applications  involving  large,  transient 
magnetic  fields  because  it  eliminates  induced  cross  currents  flowing  in  the  matrix.4 * 6  A 
matrix  of  relatively  high  resistivity  copper-nickel  alloy  (cupronickel)  is  sometimes  used 
to  impede  these  cross  currents.  Typical  cross-sectional  area  ratios  of  matrix  to  super- 
conductor are  between  1 to  1 and  2 to  1 . 


Niobium  tin  which  has  a Tc  of  18K  and  an  Hc2  of  25T  is  used  in  applica- 
tions where  fields  of  7T  (the  upper  limit  of  NbTi  for  practical  applications)  or  more 
are  required.  However,  Nb3Sn  is  very  brittle  and  cannot  be  drawn  into  wire  form. 
Until  recently,  Nb3Sn  has  been  available  only  in  rectangular-cross-section  tape,  or 
ribbon.  This  ribbon  consists  of  Nb3Sn  a few  /urn  thick  laid  on  metallic  substrate  with  a 
normal  stabilizer  as  the  outside  layer:  this  construction  allows  the  composite  to  be  bent 
without  damaging  the  superconductor.  Coils  made  from  Nb3Sn  tape  are  usually  of 
the  pancake  variety.  A pancake  is  simply  a length  of  tape  wound  on  a spool  just  as  a 

4 Z.  J.  J.  Stckly,  “State  of  the  Art  of  Superconducting  Magnets,”  / Appl.  Phys. . 42,  p.  65.  1971. 

•’  M.  N.  Wilson,  ct  al.,  “Experimental  and  Theoretical  Studies  of  Filamentary  Superconducting  Composites." 
/.  Phys.,  D,  Vol.  3,  No.  1 1,  1970. 

6 R.  R.  Critchlow,  B.  Zcitlin,  and  E.  Gregory,  “The  Effect  of  Twist  on  AC  Loss  and  Stability  in  Multistrand  Super- 
conducting Composites,”  Appl.  Phys.  Lett.  Vol.  15,  No.  7,  1969. 


r-  S4V*  t - -m 


roll  of  scotch  tape.  To  make  a solenoid,  many  of  these  pancakes  are  stacked  on  top 
of  one  another  and  connected  in  series.  Solenoids  made  from  Nb3Sn  tapes  have  pood 
high-field  performance;  but,  not  having  the  advantages  of  small  filamentary  twisted 
conductors,  their  transient  field  behavior  is  generally  not  as  good  as  that  of  NbTi. 
Recently,  Nb3Sn  has  become  commercially  available  in  multifilament  form,  but  at 
present  this  conductor  is  expensive  and  has  been  used  in  few  practical  applications. 


A vanadium-gallium  alloy,  V3Ca,  is  the  most  recent  entry  into  the  field  of 
commercial  superconductors.  V3Ga  has  even  better  high-field  performance  than 
Nb3Sn;  however,  since  it  is  relatively  new  and  expensive,  it  has  not  been  used  in  many 
applications. 


II.  THEORY  OF  AC  LOSSES  IN  SUPERCONDUCTORS 

6.  Hysteresis.  Above  H j , the  flux  penetration  into  the  superconducting 
material  causes  hysteresis  losses  under  transient  field  conditions.  The  AC  creates 
alternating  fields  in  the  superconductor  so  at  high  peak  currents  the  material  is 
constantly  being  driven  around  its  magnetization  cycle.  Methods  of  calculating  these 
hysteresis  losses  have  been  published,7'10  and  many  of  these  methods  are  widely 
accepted  today.  The  critical-state  model  due  to  Bean"  12  is  one  of  the  more  widely 
used  methods  for  calculating  these  losses.  The  major  Haw  in  the  Bean  and  London13 
critical-state  models  is  that  they  assume  that  the  critical-current  density  of  the  super- 
conductor, Jc,  is  independent  of  magnetic  field  and  that  the  full-current  density  flows 
in  regions  of  the  superconductor  where  magnetic  fields  have  penetrated.  Clearly,  the 
assumption  of  the  independence  of  Jc  on  the  field  is  not  a valid  assumption.  However, 
using  the  assumptions  of  the  critical-state  model,  the  hysteresis  losses  in  supercon- 
ductors can  be  calcualted  by  any  one  of  several  fundamentally  equivalent  techniques: 

( 1 ) The  power  dissipated  can  be  calculated  from  the  area  enclosed  within 
the  magnetic  hysteresis  loop,  i.e.,  power  loss  per  cycle  is 

P a JdV  / H • dB. 
volume  cycle 

7 R.  Itancox,  “Calculation  of  AC  Losscsin  a Type  II  Superconductor,”  Proc  II I I Vol.  1 1 3,  No.  7,  1 OOf. 

Q 

II.  London,  “AC  Losses  in  Superconductors  of  the  Second  Kind.”  Phvs  Lett,  Vol.  6,  No.  2,  1963. 

9 

(..  II.  Morgan,  ‘Theoretical  Behavior  of  Multicore  Superconducting  Wire  in  a Time-Varying  Cniforni  Magnetic 
f ield.”./.  Appl.  Phvs.,  Vol.  41,  No.  9,  1970. 

10  C.  I’.  Bean,  Phvs.  Rev.  Lett..  8:  250  (1962). 

11  .... 

I hid. 

12. 

“ C».  Bogncr  and  W.  Ilein/.el,  “Alternating  Current  Losses  in  Hard  Superconductors.”  .V.  St.  Ph  i.,  Vol.  7,  pp.  91-99. 
* ^ II.  London.  “AC  Losses  in  Superconductors  of  the  Second  Kind,”  Phvs.  Lett.,  Vol.  6,  No.  2,  1963. 
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(2)  The  loss  can  be  calculated  by  integrating  the  Poynting  Vector  over  the 
sample  surface  and  over  a complete  cycle, 

P a fdt  / l£  x 11  dA, 
cycle  surface 

where  dA  is  area  element  normal  to  surface. 

(3)  The  electric  field  E is  calculated  as  a function  of  position  and  time,  and 
the  local  power  dissipation  is  calculated  by  integrating  E • J over  a cycle.  This  method 
is  quite  useful  in  applications  of  the  Bean  Model  in  which  J is  either  zero  or  Jc. 

Other  critical-state  models  including  a field  dependence  of  Jc  have  been 
developed;  however,  these  models  are  more  complex  and  will  not  be  described  in  this 
report.  Hysteresis  losses  have  been  experimentally  shown  to  vary  linearly  with  fre- 
quency and  have  revealed  power  losses  proportional  to  the  square  or  cube  of  the 
current  and  field.14 ',8 


7.  Eddy  Currents.  Another  source  of  power  losses  in  superconducting  coils  is 
the  How  of  eddy  currents  in  the  normally  conducting  stabilizer  of  superconducting 
wires.  Copper  is  the  most  common  stabilizer  material;  and  at  superconducting 
temperatures  it  has  a conductivity  100  times  greater  than  conductivity  at  room 
temperature;  this  high  conductivity  creates  very  favorable  conditions  for  the  flow  of 
eddy  currents  when  transient  magnetic  fields  are  experienced.  These  losses  arc  readily 
calculated  by  the  application  of  Faraday’s  law  to  the  geometry  of  interest.  If  the  case 
of  a multifilament  superconducting  wire  in  a sinusoidal  AC  field  is  considered,  one  can 
obtain19  an  expression  for  the  power  loss  per  unit  volume  of  stabilizer 


O 


where  0 is  the  power,  V is  the  volume,  r is  the  maximum  radius  of  eddy  current 
loop  (half  the  distance  between  S.C.  filament),  and  p is  the  matrix  resistivity.  Carrying 


14  (i  Bogner  am!  W.  Ilcinzcl.  “Alternating  Current  Losses  in  Hard  Superconductors,"  S.  Si.  Elec.,  Vol.  7.  pp.  93-99. 

' S . II.  Morgan,  cl  al..  “Measurements  of  Lnergy  Losses  in  Pulsed  Superconducting  Magnets,"  / A ppl.  Phys.. 
Vol.  40,  No.  4,  pp.  1821-1829,  1964. 

lh  W.  r.  Beall.  Ir.  and  W.  Meycrhoff,  “AC  I nergy  Losses  Above  and  Below  H j in  Ni  and  Ni7.r.”  J.  Appl.  Phys.. 
Vol.  40.  No.  5.  pp.  2052-2059,  1969. 

17  l>.  Ilo,  M.  C.  Robinson.  V.  Srivastava.  and  R.  Stevenson,  “AC  Loss  as  a I unction  of  Current  and  I stcrnal 
Magnetic  I ield  in  Commercial  NbTi  Superconductors,”  Tech  Ppl.  AI'MI.TR-75-23. 

1 O 

’ W.  J.  C.arr,  Ir..  M S.  Walker.  I).  W.  Deis,  and  J.  II.  Murphy,  “Hysteresis  Loss  in  a Multifilament  Superconductor,' 
to  be  published. 

19 

I I . Broach  and  W.  I).  Lee.  “ The  Application  of  Loss  Models  to  Superconducting  Solenoids,’  Adr.  in  Cry.  I'ng., 
Vol.  19  (1974).  (Included  as  Appendix  A.) 
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out  the  calculation  for  an  hep  filament  structure  (Appendix  A),  one  has 
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Both  theoretical  calculation  and  experimental  data  indicate  that  eddy  current  losses 
are  proportional  to  the  square  of  the  frequency  and  the  square  ol  the  transport  current 
and  field.  In  the  region  H > Hc| . eddy-current  losses  are  usually  dominated  by 
hysteresis  losses  except  in  very  high  fields  or  high  frequencies. 


111.  METHODS  FOR  MEASURING  AC  LOSSES 


8.  Discussion  of  Methods.  There  are  several  methods  which  have  been  used  in 
recent  years  to  measure  transient  losses  in  superconductors,  and  these  can  generally 
be  classed  as  either  thermal  or  electrical  in  nature.2  0 21  The  most  popular  method 
is  a thermal  method  known  as  the  calorimetric  (boil-off)  method.  In  this  technique, 
the  sample  is  immersed  in  liquid  helium,  and  the  heat  produced  by  the  losses  results  in 
the  boil  off  of  liquid  helium  which  is  measured  by  a flowmeter.  The  measured  flow 
rate  is  then  related,  through  the  heat  of  vaporization  of  helium,  to  the  power  losses  in 
the  superconducting  sample.  The  low  sensitivity  (generally,  a few  mW)  of  this  method 
usually  dictates  a sample  of  wire  at  least  a few  meters  in  length.  The  boil-off  technique 
is  especially  suited  for  high  loss  rates  and  large  samples  since  the  lack  of  sensitivity 
causes  no  problems  in  this  case. 

There  are  several  electrical  methods  commonly  used: 

( 1 ) measurement  of  the  current  and  voltage  across  the  specimen  along  with 
the  proper  phase  relationship 

(2)  mutual-inductance  techniques 

(3)  bridge  techniques 

(4)  electronic  wattmeter  methods.22 

The  electrical  methods  are  generally  more  sensitive  than  the  thermal  methods;  some 
investigators  have  claimed  1 pW  sensitivities  which  are  required  when  small  loss  rates 
are  to  be  measured.  For  example,  the  wattmeter  method  gives  a reading  of 

^ O.  Ho,  M.  C.  Robinson,  V.  Srivastava,  and  R.  Stevenson,  “AC  Loss  as  a l unction  of  Current  and  Ixternal  Mag- 
netic f ield  in  Commercial  NbTi  Superconductors,"  Tech  Rpt.  AFMI.-TR-75-23. 

S.  L.  Wipf,  “AC  Losses  in  Superconductors."  Prnc.  1 96 S Summer  Study  on  Superconducting  Devices  and 
Accelerators.  Part  //. 

0.  Ho,  M.  C.  Robinson,  V.  Srivastava,  and  R.  Stevenson,  “AC  Loss  as  a I unction  of  Current  and  external 
Magnetic  Field  in  Commercial  NbTi  Superconductors,"  Tech  Rpt.  AFML-TR-75-23. 
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P = IV  cos  <t) 


where 


1 = sample  rms  current 

V = rms  voltage  across  sample 

<t>  = the  angle  between  V and  I. 

In  general,  it  is  probably  best  to  use  a thermal  method  if  loss  rates  are  expected  to  be 
at  least  1 watt  or  several  watts  while  electrical  methods  are  most  suitable  in  measuring 
losses  in  the  range  of  milliwatts. 

IV.  THE  AC  BRIDGE  METHOD 

9.  Maxwell  Bridge.  An  AC  null-detecting  bridge  technique  was  used  in  the 
work  described  herein.  This  method  was  chosen  mainly  because  of  the  high  sensitivity 
provided  in  this  null-detecting  method.  Many  types  of  AC  bridges  were  considered, 
but  the  one  finally  chosen  was  the  Maxwell  bridge23  (Figure  3)  which  combines  high 


sensitivity  with  ease  of  operation  since  only  variable  resistors  are  needed  thus  avoiding 
the  problems  associated  with  variable  inductors  and  capacitors.  The  balance  condi- 
tions for  this  bridge  are 

V.  A.  Brown  and  B.  P.  Ramsay,  “The  Maxwell  Bridge  at  Low  Frequencies,”  Rev.  Sci.  Inst..  Vol.  20,  No.  4. 
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L|  - (4  R;  R3 


and 

R,  - Mi 
R4 

where  L,  and  R,  represent  the  inductance  and  resistance,  respectively,  of  the  leg  of 
the  bridge  containing  the  sample  coil. 

10.  Sensitivity  and  Wagner  Grounding  Device.  An  elaborate  shielding  scheme 
which  utilized  the  Wagner  grounding  device  to  increase  the  sensitivity  and  the  signal-to- 
noise  ratio  was  used.  This  technique  has  the  effect  of  keeping  the  potential  of  the 
detector  at  the  potential  of  the  surrounding  walls,  etc,  and  thus  reduces  the  charges  on 
this  part  of  the  bridge. 


c 


Figure  4.  Maxwell  bridge  with  Wagner  grounding  device. 


The  bridge  balance  condition  Z3Z3  = Z2Z4  (Figure  4)  merely  implies  that 
the  branch  points  C,  D are  instantaneously  at  the  same  potential  and  not  necessarily 
at  0 (ground)  potential.  Thus,  current  may  still  flow  even  when  the  balance  condition 


10 


is  fulfilled  via  the  earth  capacitances  from  C and  1)  to  the  observer  and  nearby  walls, 
tables,  etc.  To  eliminate  this  effect,  two  auxiliary  branches  Z5,  Z6  capable  of 
balancing  the  impedances  Z, , Z2  are  joined  as  shown  in  the  figure  and  earthed  at  their 
common  point  E.  The  earth  capacitances  from  A and  B to  E merely  shunt  Zs  Z6.  To 
eliminate  the  capacitances  between  C,  D,  and  E,  let  the  bridge  be  balanced  as  nearly 
as  possible  to  the  minimum  deflection  on  the  detector  when  joined  between  C and  D. 
Now,  transfer  the  detector  into  the  auxiliary  network  between  C and  E balancing  by 
means  of  Z5  and  Z6 . Now,  C and  E are  nearly  at  the  same  potential,  the  latter  being 
grounded.  Reverting  to  the  original  bridge,  balance  can  now  be  more  nearly  secured, 
and  repetition  of  the  process  finally  results  in  exact  balance  with  C,  D,  and  E at  the 
same  potential,  namely,  ground.  A wiring  diagram  of  the  bridge  and  shielding  circuit 
is  shown  in  Figure  5. 

Sensitivity  tests  were  run  on  the  measuring  circuit  using  low-resistance 


copper  coils  similar  to  the  coils  to  be  tested.  Using  a PAR  Lock-in  Amplifier  as  the 
detector,  the  bridge  was  found  to  have  a sensitivity  of  about  1 x 10"s  S2.  In  actual 
operation,  an  oscilloscope,  through  Lissajous  figures,  is  used  as  the  detector  for  coarse 
balancing:  and,  then,  the  Lock-in  Amplifier  is  used  for  fine  balancing. 


1 1 . The  Meaning  of  the  Balance  Condition.  The  objective  of  these  bridge 
measurements  is  to  determine  actual  power  or  resistive  losses  in  superconducting  coils. 


R R 

Thus,  what  is  actually  represented  by  the  balance  condition,  R,  = — — —1  The  four 

R4 

legs  of  the  Maxwell  bridge  can  be  represented  by  impedance  of  the  form 

Z,  = X + iY 
Z2  = R2 
Z3  = R3 
z4  = r4 

1 + icoC4  R4 

where  legs  Z2  and  Z3  are  purely  resistive  while  Z,  is  a combination  of  resistance  and 
inductance  including  the  test  coil,  and  Z4  is  a combination  of  a resistor  and  capacitor. 
At  balance,  there  is  no  potential  across  the  detector  and  no  current  in  the  detector. 
Thus,  current  i,  is  common  to  Z,  and  Z2 , and  i2  is  common  to  Z3  and  Z4  ; so  we  have 

i 1 Z ! = i2  Z3 
i 1 Z 2 = uZ4 


% .* 


C o*  DOCTOR 
SHIELDING 


Figures.  Wiring  diagram. 


1: 


which  for  the  Maxwell  bridge  becomes 


(X  + iY) 


R4 

1 + io>C4  R4 


= R,R 


2 ,x3 


or 


R2R3  (1  + ice  C4R4)  = XR4  + iYR4 


Equating  real  parts  gives 


and 


R2R3  = xr4 


R R 

which  says  that  * — is  the  real  part  of  the  impedance  Z, . Thus,  we  are  measuring 
r4 

the  effective  resistance  of  the  test  coil  which  is  precisely  what  we  are  interested  in. 


A complete  electrical  schematic  of  the  bridge  circuit  with  detectors,  ampli- 
fiers, and  oscillator  is  shown  in  Figure  6.  The  source  of  AC  power  is  a variable  fre- 
quency sine  wave  oscillator  which  is  output  to  a differential  amplifier  where  the  signal 
is  amplified  and  isolated  without  introducing  a phase  shift.  An  oscilloscope  is  used  as 
the  detector  to  coarsely  balance  the  bridge,  and  a PAR  Lock-in  Amplifier  provides  null 
signal  detection  for  the  fine  balancing  required  for  the  10"s  S2  bridge  sensitivity.  The 
sample  coil  and  leads  are  represented  by  R,,  L,  while  the  bridge  components  R2,  R3, 
R4,  Rs,  C4,  and  C6  are  General  Radio  decade  resistors  and  capacitors.  The  leads 
connecting  the  sample  coils  with  the  remainder  of  the  bridge  circuit  were  made  of 
Evenohm  wire  so  that  the  lead  resistance  would  be  independent  of  the  liquid  helium 
level. 


V.  SAMPLES 

12.  Geometry  and  Construction.  The  wires  were  tested  in  the  form  of  small, 
multilayer  coils  wound  on  coil  forms  as  shown  in  Figure  7.  These  are  small  coils  but 
are  considerably  larger  than  traditional  small  samples  (lengths  of  wire  of  a few  centi- 
meters to  a few  meters  in  length)  with  the  typical  sample  having  wire  lengths  of  50 
meters.  All  coils  were  wound  on  identical  coil  forms  and  in  the  same  inductive  manner 
so  that  the  physical  properties  of  the  coils  were  as  similar  as  possible. 
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14 


DIMENSIONS  IN  INCHES 

Figure  7.  Sample  coil  form. 


13.  Conductor  and  Coil  Characteristics.  Tables  1 and  2 give  the  properties  of  the 
multifilament  wires  and  the  sample  coils. 


VI.  EXPERIMENTAL  METHOD  AND  DATA  REDUCTION 

14.  Taking  Data.  Exploratory  data  runs  indicated  that  the  sensitivity  of  the 
apparatus  was  such  that  lead-resistance  changes  due  to  change  in  room  temperature 
could  be  detected.  To  take  this  into  account,  the  time  was  recorded  each  time  a data 
point  was  taken.  Due  to  the  cycling  of  the  air  conditioner  in  the  lab.  the  room 
temperature  was  closely  correlated  with  the  cycle  time  of  the  air  conditioner  so  the 
time  was  used  to  make  corrections  with  the  cycle  time  of  the  air  conditioner.  This 
technique  will  be  explained  later  but  simply  amounts  to  making  linear  interpolations 
between  data  points  at  different  times. 

The  following  procedure  was  used  to  insure  that  the  resistive  and  inductive 
losses  could  be  independently  measured: 

( 1 ) Approximately  balance  bridge  on  10-  or  20-mV  scale. 
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(2)  Change  the  sensitivity  setting  on  the  Lock-in  by  two  orders  of  mag- 
nitude to  the  10-mV  or  higher  scale. 

(3)  Introduce  a large  unbalance  in  the  bridge  by  changing  R4  until  a mid- 
scale reading  is  obtained  on  the  Lock-in. 

(4)  Adjust  the  phase  control  until  a maximum  deflection  is  observed.  The 
phase  can  be  changed  by  90°  by  the  outer  ring  (quadrant  switch)  surrounding  the 
phase  shift.  When  this  is  done,  the  indicator  should  read  0,  and  when  the  quadrant 
switch  is  changed  back  to  the  0°  position  the  maximum  deflection  should  be  observed. 

1 5)  Change  R4  back  to  its  value  preceding  Step  (3). 

(6)  Switch  the  sensitivity  scale  to  the  highest  sensitivity  possible  for  an  on- 
scale  reading.  The  bridge  is  now  ready  to  be  balanced. 

As  mentioned  previously,  the  bridge  is  first  balanced  using  the  oscilloscope  as  the 
detector.  Then,  a more  sensitive  balance  is  obtained  by  using  the  Lock-in  Amplifier 
as  the  detector;  this  allows  the  bridge  to  be  balanced  within  a few  tenths  of  a milli- 
volt of  zero. 


15.  Anomalous  Resistance  Effect.  The  reduction  of  data  from  the  first  two 
samples  indicated  an  unexpected  and  initially  unexplained  behavior.  It  was  anticipated 
that  the  resistance  would  increase  continuously  with  frequency  from  0 Hz  since  the 
eddy-current  losses  increase  with  the  square  of  the  frequency.  However,  the  initial 
data  demonstrated  a behavior  as  shown  in  Figure  8 where  we  see  that  between  5 and 
350  Hz  the  resistance  decreases  with  increasing  frequency  reaching  a minimum  at 
about  350  Hz  at  which  point  it  begins  a rapid  increase  with  frequency. 

Initially,  this  behavior  was  attributed  to  frequency  response24  25  of  the 
bridge  associated  with  nonideal  components  such  as  inherent  inductance  and  capaci- 
tance ot  the  decade  resistors  and  the  capacitance  of  the  test  coil.  However,  a number 
of  complicated  calculations  were  performed  on  the  bridge  circuitry,  and  these  seemed 
to  eliminate  any  nonideal  characteristics  of  the  arms  of  the  bridge  as  the  source  of  the 
anomaly.  A separate  test  was  run  on  a typical  superconducting  test  coil  to  determine 
the  sell  capacitance  of  the  coil.  The  resonant  frequency  of  the  coil  was  determined  by 
using  the  circuit  shown  in  Figure  9 where  L, . R,.  and  (',  represent  the  inductance, 
resistance,  and  capacitance,  respectively,  associated  with  the  test  coil,  flic  resistance, 
K).  represents  the  resistance  of  the  coil  and  the  leads,  and  the  value  of  L,  has  been 

■>4 

V.  A.  Broun  and  B.  I’.  Ramsay.  “The  Maxwell  Bridge  at  Low  Ircqucneics."  Rev.  Sri.  Inst..  Vol.  20.  No.  4. 

■>S 

“ *!  Hague.  "Alternating  Current  Bridge  Methods."  as  revised  by  T.  R.  I oord  and  S.  T.  Machay,  1971.  Pitnum 
Press. 
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previously  determined.  The  resonant  frequency  was  found  to  be  1.75  x 10s  Hz.  and 
with  L,  = 3.3  x 10"3  H,  the  inherent  capacitance  of  the  coil  is 

C,  = -V-  = 2.52  x 10'10  F. 

L 

This  value  of  capacitance  is  several  orders  of  magnitude  too  small  to  account  for  the 
anomalous  behavior  in  the  0 1000  Hz  frequency  range. 

16.  Explanation  of  the  Effect.  The  following  derivation  considers  the  effect 
that  a change  in  inductance  during  each  cycle  of  AC  current  in  the  superconducting 
coil  has  on  the  balance  of  the  bridge. 

Consider  the  sinusoidal  driving  current  supplied  by  the  oscillator 

I = lo  Sin  (cot)  ( 1 ) 

where 

20 


I = peak  current. 


Figure  10  shows  the  actual  sinusoidal  variation  of  current  in  the  coil.  The  solid  curve 
of  f igure  10  represents  the  change  in  inductance  based  on  qualitative  arguments  (to  be 
given)  while  the  broken  curve  represents  a more  mathematically  manageable  situation. 
The  relative  magnitudes  of  the  curves  in  f igure  10  are  exaggerated  because  the  “bias” 
inductance,  L , is  very  large  compared  to  the  amplitude  of  the  sine  wave  variation 
since  we  are  dealing  with  very  low  currents.  In  fact,  the  currents  are  so  low  that  II  < 
II  . and  only  slight  penetration  is  assumed  to  be  caused  by  the  AC  current  in  the  coil. 

A qualitative  argument  will  now  be  presented  to  justify  the  assumption 
represented  in  Figure  10.  The  basic  assumption  is  that  some  flux  realignment  occurs  as 
the  sinusoidal  voltage  varies.  One  complete  cycle  is  broken  into  four  regions: 

( 1 ) As  the  voltage  increases  from  0 to  its  peak  value,  the  inductance  of  the 
coil  increases  (due  to  increasing  flux  penetration)  from  its  “bias”  value,  L().  to  some 
maximum  value  because  at  V , maximum  flux  penetration  occurs  and  thus  L is 
maximum. 

(2)  As  the  voltage  decreases  from  its  peak  value  back  to  zero,  the  flux 
penetration  goes  from  maximum  to  minimum,  thus  L goes  from  L to  L(). 

(3)  As  the  voltage  decreases  from  0 to  -Vmax,  the  inductance  increases 
from  L()  to  L . This  occurs  because  the  flux  penetration  does  not  depend  on  the 
sign  of  V (actually  H)  but  only  on  its  magnitude. 

(4)  For  the  same  reasons  as  in  (2).  when  the  voltage  goes  from  -Vmax  to  0. 
the  inductance  goes  from  Ln  to  Lo  . 

This  argument  explains  the  reasoning  behind  the  solid  curve  of  Figure  10, 
but  the  broken  curve  is  used  in  the  analysis  as  an  approximation  which  is  easier  to 
handle  mathematically.  Thus,  the  inductance  change  is  approximated  as  a sine  wave  of 
frequency  2 cot. 

So,  for  the  coil  inductance,  we  have 
L ~ L + L'  I 

O 

where  L'  is  the  magnitude  of  the  change  of  L.  In  the  full,  sine-wave  approximation, 
we  can  write 


L ~ Lo  + L'  sin  (2  cot). 


(2) 
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Now,  the  change  of  L.  witli  time  is 


= 2 wl.'  cos  ( 2 cot). 


from  I & M theory,  the  energy  stored  in  a coil  is 


1 = V2  LI2 


so  the  power  is 

dh  = L,  dl  + 1/2  ,2  dL  _ VI 
dt  dt  dt 

. V = L J-  + Vi  I -L  . (4) 

dt  dt 

Substituting  the  derivative  of  ( 1 ) and  (3)  into  (4)  gives 
V = LI  co  cos  cot  + co  IL'  cos  2 cot 

O 

= L 1 co  cos  cot  + L'  1 co  sin  2 cot  cos  cot  + co  1 L1  sin  cot  cos  2 cot.  (5) 

O O O O 

At  this  point,  (5)  should  be  put  into  the  response  equation  for  a lock-in  amplifier,  but 
first  we  take  the  voltage  and  divide  it  into  two  components  - one  in  phase  with  the 
current  and  one  at  7r/2  with  the  current. 

In  phase 
cot  = 7t/2 

From  (5):  V = coI0L' 

..  Z _ = -coL'  in  phase  with  1 


7t/2  out  of  phase 

cot  = o 
V = L 1 co 

O O 

Z = coLo  7r/2  out  of  phase  with  1 
0 = tt/2° 


So.  the  bridge  balance  equation  becomes 


Re  ( Z , ) = R:Rj  = R,  + R,  - coL' 
R4 


now 


K | = a oj2 

where  a is  just  a constant  of  proportionality.  Substituting  this  into  (7)  and  dividing  by 
to  gives 

R2R3  H 

p ' Ki 

co-L'=-*^ . (8) 

to 

If  the  right  side  of  this  equation  is  plotted  as  a function  of  to,  the  Y-axis  intercept 
should  give  L'.  If  this  is  a reasonable  value,  L'  < Lo  = 3 x 10~3  II  and,  in  fact.  L'  has 
a typical  value  of  10"s  II. 

This  theory  can  be  used  successfully  to  explain  the  anomalous  behavior 
described  in  the  previous  sections,  and  it  is  this  method  which  will  be  used  to  reduce 
the  resistance  data. 

17.  Data  Reduction.  I he  basic  reduction  of  the  direct-resistance  measurements 
consists  ot  compensating  for  changes  in  room  temperature  and  for  the  inductance 
changes  described  in  para.  16.  A frequency  ot  100  II/.  is  used  as  the  base  for  the  room 
temperature  compensation.  The  data-reduction  scheme  is  given  below: 


— - 


- f-  = x 

(4)  Calculate  Rj  = gjL'  = — - R0  for  each  frequency. 

. K4  J I 


R goL/ 

(5)  Calculate  — — for  each  frequency. 


to 


R,  -toL' 

(jO 

cept  and  a as  the  slope. 


(6)  Plot  ^1"  vs  to  on  cartesian  coordinates.  Then,  plot  L’  as  the  Y-inter- 
co 


(7)  Use  L'  determined  above  to  calculate  R,  which  is  the  coil  resistance 

.»  w lf=x 

R,  = 

(8)  Plot  R,  vs  f. 


RjRal  =X  R 

R4  J t 


+ coL'. 


Several  computer  programs  are  used  to  perform  this  scheme,  and  listings  of  these  pro- 
grams are  included  as  Appendices  B and  C. 


VII.  DATA  ANALYSIS  AND  RESULTS 


18.  Results.  The  method  described  in  paragraph  16  is  quite  successful  in 
accounting  for  the  seemingly  anomalous  behavior  of  the  effective  resistance  as  a func- 
tion of  frequency.  This  is  demonstrated  in  Figure  1 1 in  which  the  data  is  compared 
with  theoretical  calculations  based  on  the  method  described  in  paragraph  16.  This 
graph  is  typical  of  the  samples  and  shows  excellent  agreement  between  experimental 
results  and  theoretical  calculations. 


Computer  calculations  were  made  to  determine  the  magnetic-field  profiles 
of  the  coils;  typically,  the  maximum  field  obtained  in  the  test  coils  was  approximately 
.005  T indicating  operation  below  H ( . This  means  there  was  no  penetration  of  the 
magnetic  field  into  the  superconductor;  therefore,  hysteresis  losses  would  be  negligible, 
and  eddy  currents  in  the  stabilizer  were  expected  to  be  the  dominant  loss  mechanism. 
Some  data  was  taken  with  an  externally  applied  field;  however,  much  of  this  data  was 
unreliable  due  to  electromagnetic  instabilities  caused  by  the  mutual  inductance 
between  the  sample  and  field  coils. 

The  plot  of  Figure  12  shows  the  measured  resistance  of  each  sample  as  a 
function  of  frequency  and  shows  some  unexpected  results.  These  resistances  are  con- 
verted to  power  losses  in  Figure  13.  and  finally  the  normalized  power  losses  as  a func- 
tion of  current  are  plotted  in  Figure  14.  Since  the  dominant  loss  mechanism  is  thought 
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to  be  eddy  currents  in  the  copper  stabilizer,  the  normalization  procedure  is  to  divide 
the  power  losses  by  the  volume  of  copper  in  each  sample  coil. 

Sample  1 (Table  2)  had  considerably  lower  power  losses  than  any  of  the 
other  samples.  This  was  an  “old”  stock  of  85  filament  Airco  NbTi  wire  which  was 
untwisted;  one  possible  explanation  for  the  lower  loss  raie  is  that  the  copper  used  as 
stabilizing  material  was  not  very  pure  indicating  higher  matrix  resistivity  to  the  How  of 
eddy  currents.  Sample  2 which  is  the  same  material  as  Sample  1 but  is  twisted  shows 
loss  rates  closer  in  value  to  those  of  the  other  samples.  This  sample  also  shows  a “knee” 
in  the  curve  at  about  100  Hz.  possibly  indicating  eddy  current  shielding  due  to  the  twist 
of  2/in.  Samples  6 and  7 also  differ  only  by  the  amount  of  twist;  Sample  6 is  twisted 
5/in.  while  Sample  7 is  twisted  8/in.  The  power  losses  of  these  samples  are  almost 
identical  up  to  500  Hz  where  a knee  appears  in  the  curve  of  Sample  6 but  does  not 
appear  in  Sample  7 until  1000  Hz,  again  indicating  a shielding  effect  of  the  twist  rate. 


A sample  of  the  Airco  untwisted  conductor  was  used  for  taking  loss  data  in 
a .6  T applied,  external  field.  The  field  was  supplied  by  a 3.5  T Nb3Sn  solenoid  with  a 
1 .5-in.  ' ore.  The  results  of  this  data,  plotted  in  Figure  1 5.  indicate  considerably  higher 
losses  in  the  applied  field  as  expected  since  in  this  case  both  eddy  current  and 
hysteresis  mechanisms  contribute  to  the  losses. 

VIII.  CONCLUSIONS 

19.  Losses  Below  Hc, . The  results  of  the  study  described  in  this  report  indicate 
that  the  losses  in  superconducting  coils  in  applied  fields  below  HcJ  are  not  reduced  by 
the  introduction  of  twist  to  the  filaments.  I bis  was  demonstrated  in  measurements 
made  up  to  1000  Hz;  it  is,  however,  possible  that  at  higher  frequencies  and  higher 
fields  different  results  could  be  obtained.  Such  behavior  was  indicated  by  the  knees 
observed  in  the  twisted  samples  which  showed  that  shielding  may  occur  at  frequencies 
which  depend  on  the  twist  rate  and  other  sample  properties. 

IX.  SUGGESTIONS  FOR  FURTHER  STUDY 

20.  Losses  Above  H The  loss  mechanisms  in  the  region  Hc,  < II  < Hc2  are 
thought  to  be  considerably  different  than  those  in  applied  or  self  fields  below  Hcl 
In  the  region  above  H c| , the  magnetic  field  penetrates  the  superconductor  causing 
hysteresis  effects  which  can  become  the  dominant  loss  mechanism.  Such  a study  has 
been  undertaken  under  Work  Unit  109  as  well  as  studies  on  losses  at  60  Hz  AG  and 
variable  frequency  DC  up  to  10  Hz.  Although  there  is  an  abundance  of  transient-loss 
data  on  small,  non-inductively-wound  coils,  there  is  currently  little  data  on  multilayer, 
inductive  coils  which  operate  under  large  self-fields.  The  data  obtained  in  this  effort 
will  also  be  used  to  verify  models  developed  as  explained  in  Append.  A. 
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Multifilamentary  wires  of  Nb3Sn  have  become  commercially  available,26  26 
but  there  is  little  data  on  the  DC  or  AC  performance  of  eoils  wound  from  these 
materials.  Because  Nb3Sn  has  superior  high-field  properties  to  NbTi.  it  would  be  very 
useful  to  examine  the  performance  of  these  new  materials.  Since  most  potential 
military  applications  would  call  for  large  transient  fields,  an  experimental  study  ol  the 
transient  losses  would  make  a general  scientific  and  military  contribution  of  consider- 
able value. 


26  M.  Suenaga,  W.  B.  Sampson,  and  T.  S.  Luhman,  “Effects  of  Meat  Treatments  on  Superconducting  Critical 
Current  Densities  of  Multifilamentary  Nbs  Sn  Wires.  BNI,  I 7425. 

27  M.  Suenaga  and  W.  B.  Sampson,  “Superconducting  Properties  of  Multifilamentary  Nb3Sn  Made  by  a New 
Process,"  Appl.  Phvs.  Lett.,  Vol.  20.  No.  11,  lune  1972. 

28  “Industrial  News:  Multifilamentary  Nb3Sn  is  Now  Available  Commercially,”  Cryogenics.  February  1975. 

29  C.  F.  Oberly,  M.  C.  Ohmer,  and  11.  1..  Gegel,  “Properties  of  Multifilamcnt  Nb,Sn  in  Useful  Conductor  Configura- 
tions for  Superconducting  Rotating  Machines,”  Private  communication. 
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THE  APPLICATION  OF  LOSS  MODELS  TO 
SUPERCONDUCTING  SOLENOIDS 

J.  T.  Broach  and  W.  D.  Lee 

L .S’.  Army  Mobility  Equipment  Research  and  Development  Center 
Ft.  Betvoir.  Virginia 

INTRODUCTION 

The  use  of  inductors,  both  cryogenic  and  superconductive,  for  the  storage  of 
energy  in  pulsed  power  systems  of  various  types  has  recently  attracted  considerable 
attention.  The  storage  element  is  usually  operated  in  a pool  of  liquid  cryogen  but  for 
certain  applications,  integrated  refrigeration  using  helium  gas  is  desirable.  In  either 
case,  an  estimate  of  the  losses  generated  during  transients  is  an  essential  part  of  the 
design.  The  principal  loss  mechanisms  in  the  superconducting  case  tire  eddy  currents 
in  the  stabilizer  and  hysteresis  in  the  superconductor,  both  of  which  depend  upon  the 
magnitude  and  frequency  of  the  transient  field.  This  presentation  presents  the  results 
of  an  application  of  loss  models  to  a particular  configuration  of  multifilament 
superconductor.  The  eddy  current  loss  model  is  compared  separately  to  loss  measure- 
ments made  on  a coil  wound  from  commercial  copper  w ire.  Losses  calculated  using  a 
simple  linear  model  of  the  field  in  the  w inding  are  compared  w ith  those  obtained  from 
a more  sophisticated  calculation  using  digital  computer  techniques. 

| EDDY  CURRENT  LOSSES 

The  power  losses  due  to  eddy  currents  caused  by  a time-varying  self-field  in  a 
normally  conducting  solenoid  will  first  be  calculated.  The  magnetic  field  is  assumed 
to  be  axial  and  the  eddy  currents  are  assumed  to  flow  in  circular  paths  within  the 
w ire. 

Applying  the  integral  form  of  Faraday's  law  to  the  circular  geometry  of  an  eddy 
current  with  a maximum  radius  rm,  we  obtain  the  electric  field  driving  the  eddy 
currents 

E=\rmB  (1) 

From  Ohm's  law . the  power  density  can  be  written,  using  ( 1 ).  as 

A 

d = ('m’  4 p)62  (2) 


Working  in  spherical  coordinates,  taking  the  field  along  the  0 — 0 axis  so  that 
r„  = r sin  I). 


.1.1.  Broach  and  \\.  I).  l ie 


Integrating  to  get  the  a\eragc  loss  per  unit  volume 

<?  fo/o'/u"  [B^'lsiir  0)  4/>]i/l 

I 4 nrj  3 ,4' 

results  in 

<3  1 B’r..-  I ()/i  (5) 

Assuming  a sinusoidal  transport  current,  this  becomes 

Q 1 = U.;n)'r,;  20/<  (6| 

To  check  this  model,  measurements  of  the  effective  resistance  of  a copper  coil 
were  taken  at  various  frequencies  w ith  a Maxwell  bridge  ['  . To  obtain  the  frequence 
dependence,  the  dc  resistance  was  subtracted  from  the  measured  \alue  at  each 
frequency  The  measured  resistance  was  converted  to  power  loss  by  using  l:R,„ 
Values  of  the  power  loss  obtained  in  this  way  are  compared  in  Pig.  I with  the  results 
calculated  from  (6).  The  value  of  resistivity.  1.6  x 10  ff-m.  used  in  the  loss  calcula- 
tion was  determined. experimentally  at  4.2  k An  expression  for  H which  varies 
linearly  across  the  winding  was  used.  The  departure  of  the  measured  values  from  a 
straight  line  at  the  higher  frequencies  is  due  to  the  fact  that  the  experiment  was  run  at 
constant  voltage  across  the  bridge.  As  frequency  is  increased,  the  impedance  of  the 
coil  increases;  therefore  this  arm  draws  less  current  and  the  I R decreases  below  the 
expected  value. 

In  order  to  derive  an  expression  for  eddy  current  losses  in  the  copper  stabilizer 
of  multifilament  superconducting  wires,  the  method  used  in  the  solid  copper  wire 


Kig.  I Power  loss  vs.  frequency  for  copper 

coil  al  4.2  K f ig.  2.  "Unit  cell"  for  multitfiamcnt  conductor 
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derivation  was  adapted  to  a geometry  in  which  the  filaments  are  distributed  through- 
out the  wire  cross  section.  I he  edd>  currents  are  confined  to  the  regions  of  copper 
between  the  filaments  Considering  a cross  section  of  the  composite,  we  assume  a 
"unit  cell"  defined  by  superconducting  filaments  of  circular  cross  section  (see  Fig.  2). 
I he  composite  is  taken  to  consist  of  "hep"  cylinders  which  represent  the  filaments: 
the  spacing  ol  the  filaments  is  taken  as  a multiple  of  the  filament  radius,  allowing 
application  to  a variety  of  conductors.  The  integration  of  the  eddy  current  paths  was 
taken  in  the  regions  between  the  filaments  so  that  the  power  loss  per  unit  volume  of 
copper  was  calculated.  The  spacing  of  the  unit  cell  was  estimated  from  the  number  of 
filaments,  the  copper-to-supereonductor  ratio,  and  measurement  of  photomicro- 
graphs. Apply  ing  (3)  to  the  geometry  of  Fig.  2 and  integrating  y ields  the  loss  per  unit 


volume 


Q 

l 


2nB:  ["  rm4 

Fp  J,,  4 


tlx' 


(7) 


Referring  to  F ig  2.  rm  is  equal  to  one-half  the  distance  in  the  y direction  between  the 
outer  boundaries  of  the  unit  cell  taken  in  such  a way  as  to  exclude  the  region  occupied 
by  the  superconductor  filaments.  The  expression  obtained  for  a single  "unit  cell" 
is 


+ it  + iv  + 
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HYSTERESIS  LOSSES 

The  hysteresis  losses  in  a solenoid  wound  from  composite  superconductor  can 
be  estimated  using  a model  due  to  Bean  A long  solenoid  can  be  approximated 
as  a multilayer  scroll  wound  from  a sheet  of  superconducting  material.  The  local  power 
dissipation  per  unit  volume  derived  by  Bean  [•’]  is  given  as 


a,  — r \ 

4 («,-</,)  (4' 

This  expression  assumes  a simple  linear  variation  of  field  in  the  vv  inding  as  a function 
of  r and  assumes  no  variation  with  r (displacement  parallel  to  the  axis).  A coil  design 
based  on  NIOMAX  TC*  composed  of  101  turns  of  a single  layer  of  144-strand  braid 
having  u,  = 0.194  m.  = 0.202  m.  I = 0.792  m.  /.  = 0.3.  and  designed  for  «„ 

1 .4  T and  an  inductance  of  2 x 10  ' H with  10'  .1  stored  energy  resulted  in  calculated 
losses  of  1 .24  .1  cycle. 

For  many  solenoid  configurations,  the  simple  linear  field  profile  assumed  above 
is  a good  representation.  To  determine  the  error  which  might  be  introduced  in  design 
estimates  through  the  use  of  (9),  a more  refined  estimate  of  the  field  profile  for  the 
case  described  was  made  with  the  aid  of  a digital  computer.  The  algorithm  was  based 
on  the  superposition  of  field  contributions  due  to  a finite  number  of  current  loops 
The  field  was  found  to  be  approximately  linear  with  r near  the  midplane  with  the  field 
reversing  at  a point  near  the  outer  winding  edge.  Near  the  end  of  the  coil,  the  field 
was  proportional  to  r.  To  simplify  the  estimate  based  on  (9).  the  coil  was  divided 
into  regions  and  the  field  profile  in  a region  was  fitted  to  a simple  algebraic  expression. 
Calculated  hysteresis  losses  based  on  this  method  were  found  to  be  approximately 
5",,  lower  than  the  previous  estimate. 

* Imperial  Melal  Industries,  l td  , Birmingham.  United  Kingdom 
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The  method  described  above  has  been  compared  to  calorimetric  measurements 
of  solenoid  losses  made  at  low  frequencies  by  Jungst  [4],  Eddy  current  losses  were 
negligible  at  these  frequencies.  The  results  are  shown  in  Fig.  3.  There  is  fair  agreement 
between  the  calculated  and  experimental  results  in  this  frequency  range,  although  the 
slope  of  the  curve  dilTers  from  that  of  a curve  drawn  through  the  data  points.  The 
solenoid  upon  which  the  calculation  is  based  has  an  I D of  0.024  m,  an  OD  of  0.065  m, 
and  a length  of  0.053  m.  The  conductor  has  a diameter  of  0.0004  m and  is  composed 
of  61  35-//m  Nb  Ti  filaments  with  a twist  rate  of  four  turns/in.  The  copper-to-supcr- 
conductor  ratio  is  1.4: 1.  Central  field  is  5.2  T at  46.9  A. 

NOTATION 

(j,  inner  Minding  radius 

u:  outer  winding  radius 

H magnetic  flux  density 

sinusoidal  peak  lluxdensity 
H„  central  field  of  solenoid 

./  superconductor  filament  diameter 

h electric  field  strength 

/ rms  transport  current 

./  superconductor  critical  current  density 

I winding  length 

n scale  factor  (not  necessarily  integral)  for  filament  radius  specifying  spacing  of  composite 

Q power 

i/  power  densit) 

r radial  position 

KtU  effect i ve  resistance 

i m maximum  radius  of  eddy  current  loop 

rm  radius  of  wire 

f volume  available  for  eddy  current  paths 

(>rei‘k  letters 

/ ratio  of  superconductor  volume  to  winding  volume 

/»  resistivity 

r charge  tor  discharge)  time 

c/  angular  frequency 
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APPENDIX  B 

COMPUTER  PROGRAM  - L'  CALCULATION 


*A/(  -L  •/ 

- V K I T t TOT 

RF  TYPE  F ILF  NAMF : TP 1 

Ui  MEN SI  ON  T>.9  9)pPFF*99>pK3(99),Pa*99>pF(99)pR31<99)pRa|<99)pTF<99)p 
PI  ML (99 ) p R1 < 99)pRFT0*  99) p K 99 ) 

WRIT 

2 FORMAT  l 7FibE  WARE:  p// 

52H  bEFOPE  PROGRAM  IS  RON  HE  SURF  T HA  T THt  VALOIS 

20m  » F RIOOKKLINE  3A>  ,/ 

ASHaNL  KO  < L I IN  fc  37)  Ah  F CORRECT  FOR  THIS.  LA  1 A . ) 

1=  1 

K3 1 ( M > = 0 . 

WRITF* 1 » 6) 

6 FORMAT*  S1HRFAL)  I iv  T I ME  p R3p  R A#  An  IT  SENTINEL  K.YF  K IS  NOT  ZERO  / 
SClrtl  H I S REAL  LOOP  W ILL  EnO.  ALL  VALLES  ARE  FOP  F = TOO.) 

S REA  I:  (Op  10)T(  I )iH3(  I )jRF(  I )*K 

10  FiuRmAT*  3FR  • 5p  I I ) 

R2=70. 

REF ( I >=R2+R3* I )/R(( 1 ) 

1=1  + 1 

IF ( 1- 100) I 1 . 11 » 75 

11  IF < K ) 1 3.  5p  1 3 
M=0 

13  WRITE* 1 p 12) 

IS  FORMAT  ( 5SHREAD  IN  F R F CUENCY p T I ME  OF  ME  A SLR  EM  EN  Tp  CO  RR  E.  SPON  [.  I w G,  / 
5SHR3  AND  RApANO  L.  L MOST  PE  ZERO  OR  PROGRAM  WILL  STOP) 

M=M+  1 

15  R E A 0 ( 0 p 20  ) F ( M ) p T E ( M ) p R 3 1 ( i»| ) , R A 1 ( M ) p L 

SO  E 0 RMA T(F5p1pF5p1p2Fo.a,  II) 

IF  <L  ) 1 A6p  S5p  1 A6 

25  1 F( TE*rt)-T<  1 > ) 75p 30p 30 

30  J=S 

31  IF  C ! ECM  )-TC  J)  ) A5p  aOp  35 
35  J=J+1 

IF<  J-  I ) 31  . 75p  75 
AO  REFE=REF(J> 

OR  TG  50 
A5  J1=J-1 

REFL=REF(J1)+(RLF(J)-REF(J1))/(T<J)-T<J1>)*(TF(M)-T<J1)) 

50  WRITE* Ip  S5)F(M)pT(J)pKEFE 
55  FORMAT < 2F  6»  Op  F 8 • 5) 

REFT<M)=R2*R31<M)/Pa1(M) 

R1 00 TO  = 5»  70  777 
RFTO(M)=RFFT(M)-REFF+RIOOTO 
W*M)=S.*3. I a 1 5 9 * F * M ) 

R0=5. 68850 
R1MWL=RF TO*M)-RO 
R1  ML  < M ) = R 1 M WL/  V.  *.M  ) 

GO  TO  1 3 
1 A6  M=1 
'•RITE*  1 p 1 55) 

155  FORMA  T ( 5Fi  ,/// 

29HRESLLTS  FOR  L’  C FT FRM INA T I ON J , // ) 

1 60  WRITE(1p165) 

165  FORMAT*  A8MFREO0ENCY  An p.  FREOUENCY  R 1 ML  RF  TO  ) 

170  WRITE* 1 p 1 75) F*M)p  W*M)»  R1ML*M) p RFTO*M) 

175  FORMAT*  F6.0p6X,F10.0p7X,F9.6p3X,F7.5) 

Ms  M ♦ I 

IF  * R31  <M>  ) I 70 p 7 5p  1 70 
75  STOP 
1000  FNL 


APPENDIX  C 

COMPUTER  PROGRAM  - Ru  CALCULATION 

+9RITE  TR 

OIMENSIUN  1(  r l V9)<ki  I (9V)i 

OIMENSIUN  RMO(  lOO.KC-ET!  100) 

tRIUtl.R) 

3 FORMAT!  7HE*E  'xARE  : »// 

52b  PE  FORE  PRCCRAi-1  IS  RUN  RE  SURE  THAT  THE  VALUE  / 

aaHOF  F I GOTO  (LINE  3R  IS  CORRECT  FOR  THIS  DATA.,//) 

K 3 I (M)  = 0. 

1 = 1 

6 URITE<1»7) 

7 FORMAT  ! S1HRE.AU  IN  T IME.  R3»  l<A,  AND  SENTINEL  K ♦ Y F A IS  N0  T ZFRFI  / 
50  H THIS  RF  A 0 Lfi0P  FILL  FNU.  ALL  VALUES  ARE  EUR  F = I 00 . > 

S KEAD!0»10>T!I)»R3CI).RA!l>»K 

10  FORMAT!  3F'B*  5.  I I > 

R2=70. 

REF  < 1 >=R2+R3!  I >/RA! I ) 

1 = 1 + 1 

I F!  I- 100 ) I I . 11 » 75 

11  I F < K)  1 3.  5.  1 3 

M = 0 

13  URITE  C 1 , 1 2 ) 

1?  FORMAT!  52HREA0  Iim  F KEOUENCY. T IME  OE  M EA  SURFM  ENT  . C0RRESP0N  L IN  C.  / 
52HR3  ANU  R/.,ANU  L.  L MUST  bE  Z EKO  OR  PK0CRAM  FILL  STOP  > 

M = M*I 

15  R E A D < 0 » 20  >F!M)»TE!M)»R3I<M).Ra1!M).L 

20  FORMAT! F5* 1 . F 5. 1 »2F8* a,  i | > 

I F !L  ) 70#  25.  70 

25  IF < TE(H)-T( 1 ) ) 75. 30. 30 

30  J=2 

31  IE(TE<M)-T(J) >aa. AO. 35 

35  J=J+1 

I E ! J-  I ) 31 . 75.  7b 
4)  REFE=REFCJ> 

C-6  TO  50 
aa  J1=J-1 

AS  RE  FF  = REF ! J1 )+(REE(J)-REF(Jl >>/!T!J)-T!Jl >>*!TE!M>-T<J1 )) 

50  VRITE(1.55)F(M).T<J),REFE 

55  FORMA  T<2E6*0.FB*5> 

FtFT!M)-=  R2*R  31  ! M ) / R A1  ( M ) 

RI00T0=S.6A361 
R F TO  < M ) = R E F T < M ) • R E F E ♦ R I 00 TO 
FRITE1 l.60)F(M>.RFT0<M> 

60  FORMAT! 2HF  = » F 5*1»6A«5HRFT0=»F8*5> 

FRIT  E l 1 . 66 > 

86  E ORMA  T ! A2HkE  A U IN  KF!0*  Ir  8R0  15  NOT  ZERO  FXil  LOOP*) 

RE AU! 0. 65)KR0 
65  EORMATIIR) 

IF! KRO) 70. 1 3. 70 

70  M = 1 

9 R I T E ( 1 . 7 5 ) 

75  FORMAT! 5H  ,/// 

29HRESULTS  FOR  RO  LE T ERM I N A T I 0N : . / / 

3 7H  FRFC  UFNO Y R31  RA|  RETO  ./) 

71  VRITF(l»76)F<M>.R31<M).pA|!M>,RET0<M) 

76  FORMAT ! F7.0»5X,FS.2»3X,F7.3.3X,F7*5> 

M = M»  1 

I F<  F!M)  ) 71 , 72.  71 


7?  STOP 
FNU 
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